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亚 临界 雷诺 数 下 粗糙 圆柱 绕 流 数值 模拟 研究 
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摘 要 :基于 标准 k-e 模型 ,在 亚 临 界 雷诺 数 下 (Re =3 900) 对 不 同 表面 粗糙 度 绕 流 进行 
了 数值 模拟 研究 ,分析 了 粗糙 度 对 升力 系数 、 阻 力 系 数 、 涡 量 场 及 表面 压强 prae 影响 。 结 果 表 明 : 
随 着 表面 粗糙 度 的 增 大 ,阻力 系数 波动 的 峰值 频率 逐渐 单一 ,升力 系数 和 阻力 PHARMA 率 无 规 
律 变化 , 均 方 根 升力 系数 和 平均 阻力 系数 减 小 ; 尾 涡 排列 的 规则 性 降低 , 且 涡 泻 放 强度 减 小 ;流动 分 
离 点 向 圆柱 后 方 推移 。 在 圆柱 的 后 驻 点 ,大 粗糙 度 对 应 较 小 的 压力 系数 。 本 研究 成 果 有 助 于 圆柱 
绕 流 控制 技术 研究 。 
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Numerical simulation of flow around a circular cylinder 
with rough surface at subcritical Reynolds number 
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(1. Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanics , Wuhan University of Technology 


430063 Wuhan, China;2. Shanghai Ship Design and Research Institute ,201203 Shanghai , China; 
3. Institute of Advanced Structure Technology , Beijing Institute of Technology , 100081 Beijing , China) 


Abstract : Based on standard k-e , the two-dimensional flow around a circular cylinder with different surface 
roughness is numerically simulated at a subcritical Reynolds number of 3 900. The effects of surface rough- 
ness on the lift coefficient, drag coefficient, vorticity field and surface pressure distribution are analyzed. 
The results show that with the increase of surface roughness , the number of peak frequencies of drag coeffi- 
cient fluctuation decreases, the dominant frequencies of lift and drag coefficient changes irregularly , and 
the average drag coefficient and root-mean-squared lift coefficient decrease. The regularity of wake vortex 
arrangement and the strength of vortex shedding decreases. The flow separation point moves towards the 
rear of the cylinder as the roughness increases. At the back of the cylinder, lower pressure coefficients cor- 
respond to higher roughness. The results of this paper are helpful to the research of cylindrical flow control 
technology. 
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随 着 全 球 范围 内 海洋 资源 开采 的 需要 ,众多 海 
洋 工程 结构 物 陆续 产生 ,如 海上 风力 发 电 装置 .钻井 
隔 水 管 输 油 立 管 。 在 这 些 结构 物 中 ,圆柱 形 结构 占 
据 重要 地 位 。 在 实际 海洋 环境 中 作业 的 海洋 结构 
物 , 由 于 工艺 等 原因 会 存在 初始 表面 粗糙 度 ,而 随 着 
时 间 的 发 展 ,海洋 生物 的 附着 以 及 细小 沉积 物 的 堆 
积 会 改变 结构 物 的 表面 粗糙 度 。 粗 糙 度 的 改变 会 对 
水 动力 特性 及 流 场 状态 产生 影响 ,可 能 会 加 剧 结 
构 物 的 振动 ,使 结构 物 提前 破坏 。 因 此 ,有 必要 进行 
粗糙 度 对 结构 物流 场 的 影响 规律 研究 。 

关于 不 同 粗糙 度 下 的 圆柱 绕 流 问题 ,部 分 学 者 
展开 了 实验 研究 。Van Bearman 45 fg E Ms 9t E 
超 高 临界 雷诺 数 范围 内 对 固定 的 粗糙 圆柱 进行 了 风 
洞 实验 ,探究 了 粗糙 度 在 不 同 雷诺 数 下 对 升 阻力 系 
数 的 影响 ; Achenbach 等 四 使 用 滚 花 工 艺 在 圆柱 表 
面 形成 均匀 的 粗糙 度 , 重 点 对 粗糙 度 对 斯 特 劳 哈 尔 
ERRAT ETT T Ti s Achenbach!” 对 粗糙 圆柱 边界 
层 因 的 分 离 点 位 置 进行 了 研究 。 国 内 外 学 者 对 弹性 
支承 的 圆柱 与 其 周围 流 场 的 相互 作用 也 展开 了 一 定 
的 研究 “” ,探究 了 粗糙 度 对 流动 及 振动 响应 的 影 
WEN Buresti ^ 用 粗糙 度 不 同 的 砂 布 贴 在 圆柱 表面 
进香 风 洞 实验 ,发 现 超 临界 的 斯 特 劳 哈 尔 数 、 压 力 及 
陋 沁 系数 是 粗糙 程度 的 函数 。 对 于 粗糙 圆柱 的 研究 
大 泌 以 实验 的 方式 展开 ,主要 侧重 于 水 动力 特性 研 
究 本 但 由 于 实验 设备 的 局 限 性 ,实验 方法 只 能 采集 
到 岩 定 点 的 数据 ,不 能 得 到 整个 流 场 的 结果 。 随 着 
计 手 机 技术 的 发 展 ,数值 仿真 技术 被 越 来 越 多 地 应 
用 于 流体 力学 的 研究 。 

中 在 数值 模拟 研究 中 ,许多 不 同 的 洪流 模型 被 用 
于 处 理 满 流 场 问 题 。 马 瑞 轩 等 "采用 基于 SST in 
流 模型 的 自 适 应 尺度 模拟 方法 研究 了 串 列 双 圆柱 绕 
流 的 绕 流 特性 ,并 与 已 有 文献 中 的 数据 进行 了 详细 
地 比较 。 章 大 海 等 "采用 Spalart-Allmaras iii jfi fs 
型 对 大 范围 雷诺 数 下 非 对 称 粗 糙 单 圆柱 流 致 振动 进 
行 了 数值 模拟 。 杨 家 栋 " 和 Gao 45/9?" 均 采 用 
SST k-c 模型 对 两 向 自由 度 粗 糙 圆 柱 涡 激 振动 进行 
了 探究 ,讨论 了 粗糙 度 对 水 动力 系数 和 圆柱 振动 轨 
迹 等 响应 特性 的 影响 。Mallick 等 2 采用 大 涡 模拟 
对 粗糙 度 对 圆柱 表面 压强 分 布 的 影响 进行 了 探究 。 
Patil 4& 7 采用 高 度 非 线性 的 维 偏 微分 方程 研究 了 
表面 粗糙 度 对 细 长 运动 圆柱 绕 流 的 影响 。Yamagishi 
等 3 运用 KNG k-w 滑 流 模型 对 表面 具有 不 同 凹 楼 
的 粗糙 圆柱 进行 了 二 维 数值 模拟 研究 。 这 些 消 流 模 
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型 各 有 其 优势 和 不 足 。 除 上 述 漠 流 模型 外 ,在 Flu- 
ent 中 ,标准 k-e 模型 自从 被 提出 之 后 ,就 变 成 工程 
流 场 计算 中 主要 的 工具 。 它 具有 适用 范围 广 ` 经济、 
精度 合理 等 优点 。 

目前 对 粗糙 圆柱 的 研究 已 经 在 实验 和 数值 模拟 
方面 取得 了 一 定 的 进展 ,但 是 缺乏 对 升 阻力 系数 频 
率 特性 和 流 场 特性 的 细致 研究 。 本 研究 采用 标准 
e 模型 ,对 亚 临 界 雷诺 数 下 (Re =3 900) 光滑 和 粗糙 
圆柱 绕 流 进行 数值 模拟 ,对 实验 中 得 出 的 有 关 结 论 
进行 验证 ,并 对 粗糙 度 对 流 场 和 水 动力 系数 频率 特 
性 的 影响 展开 探究 。 


1 计算 模型 及 参数 


本 研究 的 主要 对 象 是 雷诺 数 Re =3 900 下 的 粗 
烽 圆 柱 绕 流 ,选择 数学 模型 为 标准 k-e im DER 

其 中 k-e 模型 为 双方 程 模型 ,是 在 单方 程 模型 
的 基础 上 ,引入 关于 滑 流 动能 耗 散 率 e 的 方程 。k-e 
模型 细 分 为 标准 、 重 整 化 和 可 实现 3 种 ,其 中 在 工程 
中 最 常用 的 是 标准 k-e 模型 。 漠 流动 能 耗 散 率 se 定 
义 为 


IC 
p V Ox, \ Ox, 
ii 51h EE ju, 定义 为 
E 
=pC 一 g 
M. pCs (2) 


标准 k-e 模型 的 输 运 方程 为 
(pk) Chu) _ [| 2.) k " 


ot Óx; 0x; Oy) 0x; 
G, +G, - pe -Yy * S, (3) 
0 ; ; 
(pe)  9(peu) _ 3 [fw Ps " 
ot 0X; 0x; Os) Ox; 


2 
CT (Cr + Cm6) -Cupo-*S, (4) 


E (s) 


ðu; ðu; \ ðU; 
G, =u -+ — | — , G, - fg, 
k | ) pg: Pr,ðx; 


= F b 
0x,  0x;j0x; 
p=- 2 y, =2peM},M, = Va — (6) 
p 
其 中 :6 由 平均 速度 梯度 引起 的 滑动 能 产生 ; C, 是 
由 浮力 引起 的 滑动 能 产生 ;7 为 流 场 温度 ; Yy 为 可 
压缩 溃 流 脉动 膨胀 对 总 的 耗 散 率 的 影响 ;p 为 流质 
密度 ;CC Ca 为 经 验 常数 ,在 本 研究 中 都 采用 
Fluent 中 默认 的 值 C, = 1. 4, C = 1. 922, Cm = 
0.09 ;c, 0, 1) 3 IB GERD 22 REFERIT I R 
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普 朗 特 系数 ,采用 Fluent 默认 值 o; 21.0.0, 21.3; 
Pr, 为 汕 动 普 朗 特 系 数 ,Pr, = 0. 85;g; 为 重力 加 速度 
在 i 方向 的 分 量 ;B 为 热膨胀 系数 ; M, 为 滑动 马刺 
数 ;a 为 声速 。 

如 图 1 所 示 , 本 研究 选取 圆柱 直径 为 D =0.01 m, 
流质 为 水 ,流速 为 0. 392 m/s, 雷 诺 数 为 Re =3 900。 
左 侧 边 界 为 来 流入 口 (velocity-inlet) ,与 圆心 距离 
10D, 右 侧 边 界 为 流体 出 口 (outflow) ,与 圆心 距离 为 
25D。 上 下 边界 为 滑 移 壁面 ,分 别 与 圆心 距离 10D， 
其 滑 移 速度 与 来 流速 度 相 等 ,为 0. 392 m/s。 由 此 可 
模拟 无 限 大 流 场 中 的 绕 流 运动 。 根 据 Zhao 等 (站 的 
研究 , 当 阻 塞 率 (圆柱 直径 与 流 场 宽度 的 比值 ) 小 于 
0.05 时 ,所 选取 的 计算 域 宽度 对 结果 影响 较 小 ,本 
研究 选取 阻塞 率 为 0.05 ,满足 计算 结果 不 受 计算 域 
宽 记 影响 的 要 求 。 

二 根 据 Schlinchting ” 的 离散 单元 法 思想 ,可 用 相 
同 或 不 同 的 粗糙 单元 来 离散 结构 物 表面 不 同 程度 的 
叮咬 或 凸 起 。 在 实验 设计 中 ,大 多 研究 采用 附着 了 
均 习 砂粒 的 粗糙 立 管 进行 实验 。 在 ANSYS/Fluent 
软件 中 提供 了 一 种 模拟 壁面 粗糙 度 的 模型 为 等 效 砂 
粒 模 型 ,如 图 2 所 示 ,该 模型 已 在 现 有 研究 中 多 次 使 


Mg ,其 中 K, 为 结构 的 物理 粗糙 高 度 。 


35D 


20D 


来 流入 口 “ 


滑 移 壁面 


图 1 计算 域 尺 十 
Fig.1 Size of computational domain 
Ls " 


U 


—— 


图 2 等 效 砂粒 模型 
Fig.2 Model of equivalent sand surface roughness 
在 数值 模拟 计算 中 ,将 粗糙 度 无 量 纲 化 ,用 Ks/ 
娓 量化 粗糙 度 的 大 小 。 根 据 实 际 海 详 工程 中 真实 结 


构 物 的 粗糙 度 分 布 ,选取 4 种 不 同 表面 粗糙 度 的 圆 
柱 结构 展开 研究 ,如 表 1 所 示 。 


| 
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表 1 4 种 不 同 粗糙 圆柱 的 表面 粗糙 度 


Tab.1 Surface roughness of four different rough cylinders 


组 别 1 2 3 4 


KRAHE Ks/D 1xl0-4 5x10? 1x107? 2x107? 


二 维 圆柱 绕 流 问题 中 主要 的 无 量 纲 数 有 雷诺 数 
Re 斯 特 劳 哈 尔 数 St, FHJI RZ CL、 阻 力 系数 Cs 和 
压力 系数 Cs, 其 定义 式 为 


(7) 
sp LU 为 来 流速 度 ; 为 圆柱 直径 ;x 是 流体 的 动 
JIRUE S, 为 涡 的 脱落 频率 。 
Fy F, 
Cp = 1 ,CL = 1 (8) 
zpUD PLD 
式 中 :Am 为 单位 长 度 圆 柱 在 顺 流向 受到 的 合力 ;证 
为 在 横流 向 受到 的 合力 。 
C, a =Po (9) 
I 1 » 
2P 


式 中 :p 为 当地 静 压 ;po 为 入 口 来 流 静 压 。 


2 数值 模拟 结果 


对 光滑 圆柱 绕 流 进行 数值 模拟 分 析 。 如 图 3 所 
示 , 圆 柱 上 下 两 侧 交 替 脱 落 出 旋转 方向 相反 的 涡 ,在 
圆柱 后 方形 成 规则 排列 的 双 列 线 涡 , 即 卡门 涡 街 现 
象 , 且 涡 脱落 模式 为 2S 模式 。 脱 落 涡 引 起 圆柱 表面 
压力 发 生 周期 性 变化 ,继而 形成 了 周期 性 变化 的 升 
力 系数 和 阻力 系数 。 本 研究 所 得 阻力 系数 均值 
Cy em 与 斯 特 劳 哈 尔 数 St. 与 他 人 数据 对 比如 表 2 
所 示 。 

表 2 Re =3 900 下 圆柱 绕 流 实验 及 计算 结果 总 结 


Tab.2 Experiment and calculation results of flow 


around a cylinder under Re 2 3 900 


文献 C p, mean St 
文献 [25 ] 1.00 0.210 
文献 [26 ] 1.05 0. 200 
本 研究 1.06 0.214 


由 图 4 可 看 出 , 升 阻力 系数 的 时 间 历 程 曲线 均 
为 幅 值 稳定 的 正弦 曲线 。 如 图 5 所 示 , 经 FFT 分 析 ， 
升力 系数 存在 单一 峰值 频率 ,为 8. 40 Hz; 阻 力 系数 
存在 2 个 峰值 频率 , 且 峰 值 频率 间 存 在 倍数 关系 ， 
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主峰 频率 为 16. 80 Hz ,次 峰 频 率 为 33. 60 Hz。 阻 力 
系数 的 主峰 频率 为 升力 系数 的 2 倍 ,说 明 在 圆柱 下 
游 两 侧 涡 交替 脱落 的 过 程 中 ,上 下 侧 涡 每 脱落 一 次 
即 在 一 个 运动 周期 内 , 顺 流向 的 阻力 发 生 两 次 变化 ， 
横流 向 的 升力 发 生 一 次 变化 。 
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Fig.3 Instantaneous isocontours of vorticity 
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图 4 升 阻力 系数 时 程 曲 线 
Fig.4 Time histories of lift and drag coefficients 
对 涡 量 图 和 水 动力 系数 的 分 析 表 明 , 本 研究 得 
到 的 雷诺 数 Re =3 900 下 的 结果 符合 圆柱 绕 流 的 一 
般 规 律 ,表明 本 研究 建立 的 计算 模型 以 及 采用 的 数 
值 模拟 方法 是 可 靠 有 效 的 。 
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Ci 峰值 


"D 0 15 20 25 30 3$ 40. 


到 5 升 阻力 系数 频谱 分 析 
Fig.5 Frequency analysis of lift and drag coefficients 

在 光滑 圆柱 的 基础 上 ,本 文 对 Re =3 900 ,粗糙 
度 1x10 K,/D«s2 x10 习 的 圆柱 绕 流 进行 了 数 
值 模拟 研究 。 

由 图 6 可 看 出 ,升力 系数 的 幅 值 和 阻力 系数 的 
均值 随 粗 糙 度 的 增 大 而 下 降 。 均 方 根 升力 系数 
Cims 和 平均 阻力 系数 Co wm 结果 如 表 3 所 示 。 从 表 
中 可 看 出 , 当 粗 糙 度 K/D =1 x 10 一 时 , C1, 和 
的 计算 结果 与 光滑 圆柱 的 计算 结果 十 分 接近 ， 
表明 小 粗 烟 度 对 圆柱 绕 流 特性 的 影响 较 小 。CL 只 
和 Cp ww 随 粗 糙 度 的 增 大 均 呈 下 降 的 趋势 ,与 Bear- 
man 457 关于 阻力 得 出 的 实验 结论 相符 。 

X3 不同 粗 糙 度 的 C1 ,和 Cn ,计算 结果 


Tab.3 Calculation results of C, „ms and C, mean 


with different surface roughness 


表面 粗糙 度 Ks/D CL ass CD mean 
0 0. 608 1.066 

1 x10-7* 0. 606 1.065 
5x107? 0.510 0.956 
1x107? 0.351 0.779 
2x107? 0. 102 0.486 


进一步 对 升 阻力 系数 的 频率 特性 进行 分 析 。 由 
图 6(a) (b), (o HE H: MHR K/D H 1x 
10 *,5 x10 ^8 1 x 10 ^t] , REJ RAI ERE FR] D Fe h 
线 具 有 半 周 期 性 , 即 波动 频率 不 单一 ,存在 两 个 较为 
明显 的 峰值 频率 。 升 力 系数 的 时 程 曲线 为 幅 值 稳定 
的 正弦 曲线 ,存在 单一 的 峰值 频率 ,依次 为 8.40 Hz, 
8.60 Hz 和 3.12 Hz。 

由 图 6(d) 可 看 出 : 当 粗 糙 度 K/D 92 x 107 
时 ,阻力 系数 的 时 间 历 程 曲线 具有 完全 周期 性 , 即 只 
存在 单一 的 峰值 频率 ,16. 60 Hz; 升 力 系数 的 时 间 历 
程 曲 线 保 持 为 正弦 曲线 ,峰值 频率 为 8. 40 Hz。 

整体 而 言 ,升力 系数 的 时 间 历 程 曲线 在 4 种 粗 
糙 度 下 均 为 幅 值 稳定 的 正弦 曲线 ,其 峰值 频率 因 粗 
糙 度 的 增加 呈现 不 规则 变化 。 随 着 表面 粗糙 度 的 增 
大 ,阻力 系数 波动 的 峰值 频率 逐渐 变 得 单一 ,其 时 程 
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曲线 逐渐 呈现 正弦 曲线 的 形态 , 且 升 阻力 系数 的 主 
= [G5 
1.0 12 
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导 频 率 比 保持 为 1: 2。 


0 0 a J cl a 407 
频率 /Hz 


5 10 15 20 25 30 35 40. 


(d) K/D-2 X10? 
图 6 不 同 粗糙 度 升 阻力 系数 时 程 曲 线 和 频谱 分 析 


Fig.6 Time histories and frequency analysis of lift and drag coefficients with different surface roughness 


图 7 为 不 同 粗糙 度 下 流 场 充分 发 展 后 的 涡 量 
图 。 由 图 可 知 ,与 光滑 圆柱 一 致 ,4 种 不 同 表 面 粗糙 
度 下 涡 脱 落 模 式 均 为 2 模式 , 随 着 粗糙 度 的 增 大 涡 
的 形态 略 有 变化 。 粗 糙 度 KM/D =1 x10… 时 圆柱 尾 
涡 形 态 与 光滑 圆柱 相似 ,再 次 验证 小 粗糙 度 对 圆柱 
绕 流 特性 的 影响 较 小 。 随 着 粗糙 度 的 增 大 , 尾 涡 形 
状 逐 渐 被 拉 长 , 尾 涡 排列 的 规则 性 降低 , 且 涡 脱落 的 
强度 减 小 。 结 合 图 6 可 知 , 涡 脱落 强度 降低 引起 升 
力 系数 幅 值 和 阻力 系数 均值 减 小 。 


4 种 粗糙 度 下 圆柱 的 平均 压力 系数 分 布 情 况 如 
图 8 所 示 。 从 图 中 可 以 看 出 ,4 种 粗糙 度 下 的 曲线 
形状 大 致 相同 。 随 着 粗糙 度 的 增 大 ,曲线 最 低 点 对 
应 的 Cos 值 减 小 , 且 该 值 对 应 的 角度 增 大 。 通 过 
对 角度 -压力 系数 曲线 进行 二 阶 求 导 得 出 粗糙 度 
Ks[D 为 1x10…、5x10 一 时 流动 分 离 点 均 为 0 = 
100? ,粗糙 度 K,/D 为 1 x10 7,2 x10 一 时 流动 分 离 
点 分 别 为 105" 110"。 由 此 可 知 , 随 着 粗糙 度 的 增 
大 ,流动 分 离 点 向 圆柱 后 方 推移 。 在 曲线 到 达 最 低 
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点 前 ,粗糙 度 对 压力 的 影响 不 大 。 在 圆柱 的 后 驻 点 ， 
大 粗糙 度 对 应 较 大 的 C, mean CABOS HEB) 。 圆 柱 
所 受阻 力 由 压 差 阻力 与 摩擦 阻力 构成 ,在 本 研究 模 
拟 的 情况 下 , 压 差 阻力 占 主导 地 位 ,在 图 8 中 , 随 着 
粗糙 度 的 增 大 ,圆柱 前 后 驻 点 压 差 减 小 ,这 与 图 6 中 
阻力 系数 减 小 的 现象 一 致 。 


(a) K/D-A X10? 


(b) K/D-5 x 10^ 


(c) K/D-1xX 10? 


(d) K/D-2 X10? 


图 7 不 同 粗糙 度 涡 量 图 


Fig.7 Instantaneous isocontours of vorticity 


with different surface roughness 
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8 不 同 粗糙 度 的 平均 压力 系数 分 布 


Fig.8 Distribution of mean pressure coefficient 


with different surface roughness 
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本 研究 基于 RANS 模拟 方法 ,对 亚 临 界 雷 诺 数 
下 光滑 和 粗糙 圆柱 绕 流 进行 了 研究 ,分析 了 不 同 粗 
糙 度 下 的 尾 流 特征 ,考察 了 粗糙 度 对 升 阻力 系数 、 涡 
量 场 及 表面 压强 分 布 情况 的 影响 ,主要 结论 如 下 。 

1) 升 阻力 系数 的 时 程 曲线 在 圆柱 光滑 和 粗糙 的 
情况 下 均 为 近似 的 正弦 曲线 。 随 着 表面 粗糙 度 的 增 
大 ,升力 系数 和 阻力 系数 的 主导 频率 呈现 不 规则 变 
化 , 均 方 根 升力 系数 和 平均 阻力 系数 减 小 。 

2) 随 着 粗糙 度 的 增 大 , 尾 涡 形 状 逐 渐 被 拉 长 , 尾 
涡 排 列 的 规则 性 降低 , 且 涡 汽 放 强度 减 小 。 

3) 随 着 粗糙 度 的 增 大 ,流动 分 离 点 向 圆柱 后 方 
推移 。 在 曲线 到 达 最 低 点 前 ,粗糙 度 对 压力 的 影响 
不 大 。 在 圆柱 的 后 驻 点 ,大 粗糙 度 对 应 较 大 的 
C, wm 值 (绝对 值 较 小 ) 。 


参考 文献 : 


[1] GAO Y,ZONG Z,ZOU L,et al. Numerical simulation of vortex-in- 


duced vibration of a circular cylinder with different surface rough- 
nesses[ J ]. Marine structures ,2018 ,57 :165-179. 

[2] VANHN P. The Reynolds number dependency of the steady and 
unsteady loading on a slightly rough circular cylinder; From subcrit- 
ical up to high transcritical flow state [J]. Journal of fluids and 
structures ,2015 ,55 :526-539. 

[3] BEARMAN P W,HARVEY J K. Control of circular cylinder flow 

by the use of dimples[ J]. AIAA journal, 1993,31 (10) : 1753- 

1756. 

[4] ACHENBACH E, HEINECKE E. On vortex shedding from smooth 

and rough cylinders in the range of Reynolds numbers 6 x 10? to 

5 x 109 [ J]. Journal of fluid mechanics ,1981 ,109 239-251. 

[5] ACHENBACH E. Influence of surface roughness on the cross-flow 

around a circular cylinder[ J ] . Journal of fluid mechanics ,1971 ,46 

(2) :321-335. 

[6] KIU K Y,STAPPENBELT B,THIAGARAJAN K P. Effects of uni- 


form surface roughness on vortex-induced vibration of towed vertical 
cylinders[ J ]. Journal of sound and vibration, 2011, 330 (20) : 
4753-4763. 

[7] ZHOUB,WANG X,GUO W ,et al. Experimental study on flow past 
a circular cylinder with rough surface [ J]. Ocean engineering, 
2015,109 7-13. 

[8] OKAJIMA A,NAGAMORI T, MATSUNAGA F,et al. Some experi- 
ments on flow-induced vibration of a circular cylinder with surface 
roughness [ J ]. Journal of fluids and structures, 1999, 13 (7) : 
853- 864. 

[9] GAO Y,LIU L M,ZOU L,et al. Effect of surface roughness on vor- 


tex-induced vibration of a freely vibrating cylinder near a stationary 


微 信 公众 号 :应 用 力学 学 报 


980 应 


plan wall[ J]. Ocean engineering ,2020 ,198 :106837. 

[10] GAO Y,FU S, WANG J,et al. Experimental study of the effects of 

surface roughness on the vortex-induced vibration response of a 

flexible cylinder[ J]. Ocean engineering. 2015 ,103 :40-54. 

[11] BERNITSAS M M,RAGHAVAN K. Reduction/suppression of VIV 

of circular cylinders through roughness distribution at 8 x 10° < Re < 

1.5 x 10% [ C ]//Proceedings of ASME 2008 27th International 

Conference on Offshore Mechanics and Arctic Engineering. [ S. 1. ] : 

American Society of Mechanical Engineers ,2008 : 1001-1005. 

[12] BERNITSAS M M, RAGHAVAN K,DUCHENE G. Induced separa- 

tion and vorticity using roughness in VIV of circular cylinders at 

8 x10? < Re «1. 5 x 1 [ C ]//Proceedings of ASME 2008 27th 

International Conference on Offshore Mechanics and Arctic Engi- 

neering. [ S. 1. | : American Society of Mechanical Engineers ,2008 : 

993-999. 

[13] ZHOU B,WANG X,GUO W ,et al. Control of flow past a dimpled 
s circular cylinder [ J ]. Experimental thermal and fluid science, 

52015 ,69(6) 19-26. 

[È= zuoU B, WANG X,GUO W,et al. Force and flow characteristics 
of a circular cylinder with uniform surface roughness at subcritical 
Reynolds numbers| J]. Applied ocean research ,2015 ,49 :20-26. 

GB munem G. The effect of surface roughness on the flow regime a- 

"round circular cylinder[ J ] . Journal of wind engineering and indus- 
trial aerodynamics,1981,8 (1) :105-114. 

[1611 马 瑞 轩 , 李 文 建 , 左 筷 成 . 基于 SST-SAS 方法 与 FW-H 方程 的 
EN 列 双 圆柱 数值 模拟 研究 [J]. 应 用 力学 学 报 . 2020,37 (3) : 
e 1057-1064. 

.二 MA Ruixuan,LI wenjian, ZUO kongcheng. Numerical simulation of 

5 flows around tandem cylinder based on SST-SAS method and FW-H 

< equation[ J]. Chinese journal of applied mechanics ,2020 ,37 (3 ) : 
1057-1064 ( in Chinese) . 

[VIE SESS , EGO ERY , SE. 非 对 称 粗糙 带 对 单 圆柱 流 致 振动 
T 特性 影响 研究 [中 .船舶 力学 ,2019 ,23 (10) :1177-1186. 


© 


过 


投稿 网 站 :http ;// cjam. xjtu. edu. en 


用 力学 学 报 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


ChinaXiv& fERBTI a 39 æ 


ZHANG Dahai , WANG Wenhao, LI Tianjiao, et al. Research on flu- 
id-induced vibration of single circular cylinder with asymmetrical 
roughness strip[ J ]. Journal of ship mechanics, 2019, 23 ( 10) : 
1177-1186 (in Chinese). 


] 杨 家 栋 . 粗糙 圆柱 体 涡 激 振动 响应 数值 模拟 [D] . 成 都 :西南 石 


油 大 学 ,2017. 

GAO Y,LIU L,ZOU L,et al. Effect of surface roughness on vortex- 
induced vibrations of a freely vibrating cylinder near a stationary 
plane wall[ J]. Ocean engineering,2020 ,198 ; 106837. 

GAO Y,ZHANG Z, ZOU L, et al. Effect of surface roughness and 
initial gap on the vortex-induced vibrations of a freely vibrating cyl- 
inder in the vicinity of a plane wall[ J]. Marine structures, 2020 , 
69 : 102663. 

MALLICK M, KUMAR A, MURMU A. Flow modeling in various 
cylindrical surfaces[ J]. Aquatic procedia ,2015 ,4 :834- 840. 
PATIL P M, HIREMATH P S. Effects of surface roughness on 
mixed convection nanoliquid flow over slender cylinder with liquid 
hydrogen diffusion [ J ]. International journal of hydrogen energy, 
2019 ,44(21) :11121-11133. 

YAMAGISHI Y , OKI M. Numerical simulation of flow around a cir- 
cular cylinder with curved sectional grooves[ J ]. Journal of visual- 
ization , 2007 ,10(2) :179-186. 

ZHAO M, KAJA K, XIANG Y, et al. Vortex-induced vibration 
(VIV) of a circular cylinder in combined steady and oscillatory 
flow[ J]. Ocean engineering ,2013 (8) :7383-7395. 

SCHLICHTING H. Boundary-layer theory [ M ]. New York ; Mcgraw- 
hill Book Company ,1979. 


FRANZINI G R, GONCALVES R T, MENEGHINI J R, et al. Ex- 
perimental investigation into the flow around a stationary and yawed 
cylinder under asymmetrical end conditions [ J]. International jour- 
nal of offshore & polar engineering,2014 ,24 (2) :90-97. 

(编辑 S 


微 信 公众 号 :应 用 力学 学 报 


